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applied to the data before directly recycling. But it is important 
that such an optical, classical finite-state machine is demon- 
strated to prove the feasibility of parallel-optical digital logic. 


Conclusions 


The operation of a number of digital all-optical loop circuits, 
based on bistable nonlinear interference filter (NLIF) devices, 
has proven the capability of all-optical switches of this type to 
form systems of extensible optical restoring logic. Furthermore 
the successful switching of cascaded arrays of such all-optical 
digital gates has confirmed their anticipated potential for creat- 
ing parallel processing systems. 

We have chosen to use electronic methods for control and 
clocking of these circuits. There is nothing ‘optically impure’ in 
this approach in the régime where switch times are greater than 
total-loop time. Indeed, it provides an important interface with 
conventional electronics. For various logic decision purposes 
one could anticipate that the ‘optical classical finite-state- 
machine’, as depicted, could be added to an electronic machine 
to multiply the processor power by the number of parallel optical 
channels (single instruction multiple data-stream). 

The interference film devices that have been used lend them- 
selves to parallelism because they are intrinsically uniform to a 
high degree. The pass-band (~4 nm FWHW) has been measured 
to be constant in width and centre wavelength to <3% over an 
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area of ~25 mm? with a spatially mapping Fourier transform 
interferometric spectrometer. 

The effects of thermal diffusion on spot-size scaling, gate 
cross-talk and switching speeds have been studied both theoreti- 
cally and experimentally’''*. Sample pixellation is essential to 
prevent cross-talk in large arrays. We are pursuing fabrication 
techniques that include pixellation of the active layers and of 
thermally designed substrates. Even at the present switching 
power levels, 1 mW for 10 1m spot sizes and 30 ps switching 
times, the achievable parallelism is 10° gates per watt. For 
pixellated thermal devices optimized for low-energy switching 
the predicted parameters are 1 pJ ym~? (for example, 100 W, 
10 pm’, 100 ns)'*?!, This would yield up to 10'! Hz W_'. Given 
easily achievable heat sinking of 10 W cm’ this implies a switch- 
ing capacity of 10’? Hz cm~2. Such performance would permit 
experimental assessment of a wide variety of all-optical parallel 
information processing in the near future. 

We acknowledge the essential contributions of J. G. H. 
Mathew, M. R. Taghizadeh, N. Craft, I. Redmond and R. J. 
Campbell to the experimental work. Partially supported by the 
UK Science and Engineering Research Council, Joint Opto- 
Electronic Research Scheme (JOERS); the work was carried 
out within the framework of the European Joint Optical Bistabil- 
ity Project (EJOB) of the Commission of the European Com- 
munities. 


13. Wherrett, B. S., Hutchings, D. & Russell, D. J. opt. Soc. Am. B 3, 351-362 (1986). 

14. Janossy, I., Mathew, J. G. H., Abraham, E., Taghizadeh, M. R. & Smith, S. D. IEEE J. 
Quantum Electron. (in the press). 

15. Mathew, J. G. H., Taghizadeh, M. R., Abraham, E., Janossy, I. & Smith, S. D. Optical 
Bistability LH (eds Gibbs, H. M., Smith, 8. D., Mandel, P. & Peyghambarian, N.) 57-60 
(Springer, Berlin, 1986). 

16. Janossy, I., Taghizadeh, M. R., Mathew, J. G. H. & Smith, S. D. IEEE J. Quantum Electron. 
21, 1447-1452 (1985). 

17. Walker, A. C. Optics Commun. 59, 145-150 (1986). 

18. Gibbs, H. M., Tarng, S. S., Jewell, J. L., Weinberger, D. A. & Tai, K. Appl Phys. Lett, 41, 
221-223 (1982). 

19. Ovadia, S., Gibbs, H. M., Jewell, J. L., Sarid, D. & Peyghambarian, N. Opt. Engy 24, 
565-568 (1985). 

20. Miller, D. A. B., Henry, J. E., Gossard, A. C. & English, J. H. Appl Phys. Lett. 49, 821-823 
(1986). 

21. Wherrett, B. S. in From Optical Bistability Towards Optical Computing (Proc. European 
Joint Optical Bistability Project) (eds Smith, S. D., Mandel, P. & Wherrett, B. S.) Ch. 
5.7 (North Holiand, Amsterdam, in the press). 








Mitochondrial DNA and human evolution 
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Mitochondrial DNAs from 147 people, drawn from five geographic populations have been analysed by restriction mapping. 
All these mitochondrial DNAs stem from one woman who is postulated to have lived about 200,000 years ago, probably 
in Africa. All the populations examined except the African population have multiple origins, implying that each area was 


colonised repeatedly. 


MOLECULAR biology is now a major source of quantitative and 
objective information about the evolutionary history of the 
human species. It has provided new insights into our genetic 
divergence from apes!-® and into the way in which humans are 
related to one another genetically*'*. Our picture of genetic 
evolution within the human species is clouded, however, because 
it is based mainly on comparisons of genes in the nucleus. 
Mutations accumulate slowly in nuclear genes. In addition, 
nuclear genes are inherited from both parents and mix in every 
generation. This mixing obscures the history of individuals and 


“Present address: Department of Genetics, University of Hawaii, Honolulu, Hawaii 96822. 


allows recombination to occur. Recombination makes it hard 
to trace the history of particular segments of DNA unless tightly 
linked sites within them are considered. 

Our world-wide survey of mitochondrial DNA (mtDNA) adds 
to knowledge of the history of the human gene pool in three 
ways. First, mtDNA gives a magnified view of the. diversity 
present in the human gene pool, because mutations accumulate 
in this DNA several times faster than in the nucleus'*. Second, 
because mtDNA is inherited maternally and does not recom- 
bine”, it is a tool for relating individuals to one another. Third, 
there are about 10'° mtDNA molecules within a typical human 
and they are usually identical to one another!”~*. Typical mam- 
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Fig. 1 Histogram showing the number of site differences between 
restriction maps of mtDNA for all possible pairs of 147 human 
beings. 


malian females consequently behave as haploids, owing to a 
bottleneck in the genetically effective size of the population of 
mtDNA molecules within each oocyte’. This maternal and 
haploid inheritance means that mtDNA is more sensitive than 
nuclear DNA to severe reductions in the number of individuals 
in a population of organisms'>. A pair of breeding individuals 
can transmit only one type of mtDNA but carry four haploid 
sets of nuclear genes, all of which are transmissible to offspring. 
The fast evolution and peculiar mode of inheritance of mt DNA 
provide new perspectives on how, where and when the human 
gene pool arose and grew. 


Restriction maps 


MtDNA was highly purified from 145 placentas and two cell 
lines, HeLa and GM 3043, derived from a Black American and 
an aboriginal South African (!Kung), respectively. Most placen- 
tas (98) were obtained from US hospitals, the remainder coming 
from Australia and New Guinea. In the sample, there were 
representatives of 5 geographic regions: 20 Africans (represent- 
ing the sub-Saharan region), 34 Asians (originating from China, 
Vietnam, Laos, the Philippines, Indonesia and Tonga), 46 
Caucasians (originating from Europe, North Africa, and the 
Middle East), 21 aboriginal Australians, and 26 aboriginal New 
Guineans. Only two of the 20 Africans in our sample, those 
bearing mtDNA types 1 and 81 (see below) were born in sub- 
Saharan Africa. The other 18 people in this sample are Black 
Americans, who bear many non-African nuclear genes probably 
contributed mainly by Caucasian males. Those males would not 
be expected to have introduced any mtDNA to the Black 
American population. Consistent with our view that most of 


Table 1 MtDNA divergence within and between 5 human populations 


% sequence divergence 


Population 1 2 3 4 5 
1. African 0.47 0.04 0.04 0.05 0.06 
2. Asian 0.45 0.35 0.01 0.02 0.04 
3. Australian 0.40 0.31 0.25 0.03 0.04 
4. Caucasian 0.40 0.31 0.27 0.23 0.05 
5. New Guinean 0.42 0.34 0.29 0.29 0.25 





The divergence is calculated by a published method”*. Values of the 
mean pairwise divergence between individuals within populations (6,) 
appear on the diagonal. Values below the diagonal (5,,)) are the mean 
pairwise divergences between individuals belonging to two different 
populations, X and Y. Values above the diagonal (6) are interpopulation 
divergences, corrected for variation within those populations with 
equation (1). 
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Fig. 2 Sequence divergence within 5 geographic areas for each 
of 7 functional regions in human mtDNA. Sequence divergence 
(8,) was estimated from comparisons of restriction maps”®. Sym- 
bols for the 5 races are: @, Africa; O, Asia; A, Australia; QO, 
Caucasian; A, New Guinea. Along the horizontal axis are the 
numbers of restriction sites in each functional region (D loop, 
transfer RNA genes, 128 and 16S ribosomal RNA genes, NADH 
dehydrogenase subunits 1-5, cytochrome oxidase subunits and 
other protein-coding regions). 


these 18 people are a reliable source of African mtDNA, we 
found that 12 of them bear restriction site markers known”’ to 
occur exclusively or predominantly in native sub-Saharan 
Africans (but not in Europeans, Asians or American Indians 
nor, indeed, in all such Africans). The mtDNA types in these 
12 people are 2-7, 37-41 and 82 (see below). Methods used to 
purify mtDNA and more detailed ethnographic information on 
the first four groups are as described'”*; the New Guineans 
are mainly from the Eastern Highlands of Papua New Guinea”. 

Each purified mtDNA was subjected to high resolution map- 
ping’?-** with 12 restriction enzymes (Hpal, Avall, FnuDII, 
Hhal, Hpall, Mbol, TagI, Rsal, HinfI, Haelll, AluI and 
Ddel). Restriction sites were mapped by comparing observed 
fragment patterns to those expected from the known human 
mtDNA sequence”’. In this way, we identified 467 independent 
sites, of which 195 were polymorphic (that is, absent in at least 
one individual). An average of 370 restriction sites per individual 
were surveyed, representing about 9% of the 16,569 base-pair 
human mtDNA genome. 


Map comparisons 


The 147 mtDNAs mapped were divisible into 133 distinct types. 
Seven of these types were found in more than one individual; 


Table 2 Clusters of mtDNA types that are specific to one geographic 





region 
Number of Mean pairwise Average 
Geographic region-specific distance within age of 
region clusters clusters* clusterst 
Africa 1¥ 0.36 90-180 
Asia 27 0.21 53-105 
Australia 15 0.17 43-85 
Europe 36 0.09 23-45 
New Guinea 7 0.11 28-55 


* For clusters represented by two or more individuals (and calculated 
for individuals, not for mtDNA types) in Fig. 3. 

+ Average age in thousands of years based on the assumption that 
mtDNA divergence occurs at the rate of 2-4% per million years'**°. 

+ Assuming that Africa is the source, there is only one African cluster. 
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no individual contained more than one type. None of the seven 
shared types occurred in more than one of the five geographic 
regions. One type, for example, was found in two Australians. 
Among Caucasians, another type occurred three times and two 
more types occured twice. In New Guinea, two’ additional types 
were found three times and the seventh case involved a type 
found in six individuals. 


A histogram showing the number of restriction site differences ~ 


between pairs of individuals is given in Fig. 1; the average 
number of differences observed between any two humans is 9. 5, 
The distribution is approximately normal, with an excess of 
pairwise comparisons involving large numbers-of differences. 
. From the number of restriction site differences, we estimated 
the extent of nucleotide sequence divergence”® for each pair of 
individuals. These estimates ranged from zero to 1. 3 substitu- 
tions per 100 base pairs, with an average sequence divergence 
of 0.32%, which agrees with that of Brown’, who examined 
only 21 humans. 

Table 1 gives three measures of sequence divergence within 
and between each of the five populations examined. These 
measures are related to one another by equation (1): 


5 = 6,,—0.5 (6, +5y) (1) 


where 6, is the mean pairwise divergence (in percent) between 
individuals within a single population ), 6, is the correspond- 
ing value for another population (Y), 6 xy is "the mean pairwise 
divergence between individuals belonging to two different popu- 
lations (X and Y), and 6 is a measure of the interpopulation 
divergence corrected for intrapopulation divergence. Africans 
as a group are more variable (6, = 0.47) than other groups. 
Indeed, the variation within the African population i is as great 
as that between Africans and any other group ( By = =0.40-0.45). 
The within-group variation of Asians (6, = 0.35) is also compar- 
able to that which exists between groups. For Austfalians, 
Caucasians, and: New Guineans, who. show ‘nearly identical 
amounts of within-group variation (6, = 0.23-0.25), the variation 
between groups slightly exceeds that within groups. 

When the interpopulational distances (8,y) are corrected for 
intrapopulation variation (Table 1), they. become very small 
(6 = 0.01-0.06). The mean value of the corrected distance among 
populations (6 =0.04) is less than one-seventh of the mean 
distance between individuals within ‘a population (0.30).. Most 
of the mtDNA variation in the human species is therefore shared 
between populations. A more detailed analysis. supports this 
view”’ 


Functional constraints 


Figure 2 shows the sequence divergence (6,) calculated for each 
population across seven functionally distinct ee of the 
mtDNA genome. As has been found before’*?’®, the most 
variable region is the displacement loop (6,=1. 3), the major 
noncoding portion of the mtDNA molecule, and the least.vari- 
able region is the 16S ribosomal RNA gene (dx: =0. 2). In general, 
Africans are the most diverse and Asians the next most, across 
all-functional regions. 


Evolutionary tree 


Atree relating the 133 types of human mtDNA and the reference 
sequence (Fig. 3) was built by the parsimony method. To inter- 
pret this tree, we make two assumptions, both of which have 
extensive empirical support: (1) a strictly maternal.mode of 
mtDNA transmission (so that any variant’ appearing ina group 
of lineages must be due to a mutation occurring in the ancestral 
lineage and not recombination between maternal and paternal 
genomes) and (2) each individual is homogeneous for its multi- 
ple mtDNA genomes. We can ‘therefore view. the tree as a 
genealogy linking maternal lineages i in modern human popula- 
tions to a common ancestral female (bearing mtDNA type a). 
- Many trees of minimal or néar-minimal length can be made 
from the data; all trees that we have have examined share the 
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following features with Fig. 3. (1) two primary branches, one 
composed entirely of Africans, the other including all 5 of the 
populations studied; and (2).each population stems from multi- 
ple lineages connected to the tree at widely dispersed positions. 
Since submission of this manuscript, Horai et al built a tree 
for our samples of African and Caucasian populations and their 
sample of a Japanese population by another method; their tree 
shares these two features. 

Among the trees investigated was one consisting of five 
primary branches with each branch leading exclusively to one 
of the five populations. This tree, which we call the population- 
specific tree, requires 51 more point mutations ‘than does the 
tree of minimum length in Fig. 3. The minimum-length tree 
requires fewer changes at 22 of the 93 phylogenetically-informa- 
tive restriction sites than does the population-specific tree, while 
the latter tree required fewer changes at four sites; both trees 
require the same number of changes at the remaining 67 sites. 
The minimum- length tree is thus favoured by a score of 22 to 
4. The hypothesis that the two trees’ are equally compatible | with 
the data ‘is’ statistically rejected, since 22:4-is significantly 
different from the expected 13: 13. The minimum- length tree is 
thus significantly more parsimonious than the population- 
specific tree. 


African origin 


We infer from the tree of minimum length (Fig. 3) that Africa 
is a likely source of the human mitochondrial gene pool. This 
inference comes from the observation that one of the two primary 
branches leads ‘exclusively to African mtDNAs (types 1-7, 
Fig. 3) while the second primary branch also leads to African 
mtDNAs (types 37-41, 45, 46, 70, 72, 81, 82, 111 and 113). By 
postulating that the common ancestral mtDNA (type in Fig. 3) 
was African, we minimize the number of intercontinental migra- 
tions needed. to account for the géographic distribution of 
mtDNA types. It follows ‘that b i is‘a likely common ancestor of 
all non- -African and many ‘African mtDNAs (types 8-134 in 
Fig. 3). . 


Multiple lineages per race 


The second implication of the tree (Fig. 3) —that each non- 
African population has multiple origins—can be illustrated most 
simply with the New Guineans. Take, as an example, mtDNA 
type 49, a lineage whose nearest relative is not in New Guinea, 
but in Asia (type 50). Asian lineage 50 is closer genealogically 
to this New Guined lineage than to other Asian mtDNA lineages. 
Six other lineages lead exclusively to New Guinean mtDNAs, 
each originating at a different place in the tree (types 12, 13, 
26-29, 65, 95 and 127-134 in Fig. 3). This small region of New 
Guinea (mainly the Eastern Highlands Province) thus seems to 
havé been colonised by at Teast seven maternal lineages (Tables 
2 and 3).' 

In the same way, we calculate the minimum numbers of female 
lineages that colonised ‘Australia, Asja and Europe (Tables 2 
and 3). Each estimate is’ based on the number of region-specific 
clusters in the tree (Fig.3, Tables 2 and 3). These numbers, 
ranging from 15 to 36 (Tables 2 and 3), will probably rise as 
more types. of human mtDNA a are discovered. ° 


Tentative time scale 


A time scale can be affixed to the tree in Fig. 3 by assuming that 
mtDNA sequence divergence accumulates at a constant rate in 
humans. One way of estimating this rate is to consider the extent 
of differentiation: within Clusters specific to New Guinea (Table 
2; see.also refs 23 and 30), Australia®® and the New World”). 
People colonised these regions Telatively recently: a minimum 
of 30,000 years ago for New Guinea™’, 40,000 years | 280 for 
Australia®?, and 12;000 years ago for the ‘New World™.. These 
times enable us to. calculate that the mean rate. of mtDNA 
divergence within humans lies between two and. four percent 
per million years; ‘a detailed account of this calculation appears 
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D LOOP: 8 2 5 8 9 23 81-83 85 86 118; 641 (16494) 26 ; 1341 3; 207h 128-134; 
255t 39; 259a 106; 340] 112 12S rRNA: 663e 43 48,53; 712) 90; 740] 16-19; 
748 49;1240a 7-9 37 110; 1403a (1448) 11-18 20-29 34 49 50 55-57; 1463e 14; 
14848 2 4 5 811 14-17 19 23-25 30 32-35 41-43 45 52-58 61 107 110 120; 
1536f 31 VAL tRNA: 1610a 2 3 44 45 104; 1637c 47 102 108; 1667c 47 16S 
rRNA: 1715c 22 23 27 123;1917a 70; 2208a 87; 2223a (2635) 37; 2384a (2472) 
106; 2390] 2 3 8 9 49 79 93; 2734a 18 22 54 73 101; 2758k 1-8; 2849k 112; 
3123k 119 ND 1: 3315e 11;3337k 94 118 133; 3391e 41 105 122; 3537a 44; 
3592n 2-7 9 37-41; 3698f 130-132; 3842e 44; 3849e 5 46; 38991 14 18 22 23 32 
71 92 102 110 120; 3930c 110; 39441 12; 4092g 112 METIRNA: 4411a 23 35 
36-57; 4464k 129 ND 2: 4487b (70933) 111 112; 4631a 2; 4643k 87 100; 4732k 
115; 4769a 110; 4793e 76 85; 5176a 21 25 54-57; 5261e 118; 52691 71; 5357f 
17-19 TRP tRNA: 5538 130; 5552c 10 L. ORIGIN: 5742] 15 44; 5754i ($755) 
81 CO 1: 5978a 52; 5984b/5983g 87; 5985k/5983g 133 134; 5996a 87; 6022a 74; 
6166f (6168) 63; 62119 37-41 82; 6260e 95; 6356c 5 9; 6377c 68; 64091 (7854) 
68; 65071 78; 66709 12; 6699b (8719,8723) 78; 68719 18 19; 6915k 5 49; 6931g 
48; 6957e 44 45; 7025a 75-86 93 94 103-112; 7055a 3 6 7; 7241k 72; 73351 6 7 
30 92 120; 7347e 124 SER tRNA: 74611 71; 7474a 121 COll: 7677 55 73; 
7750c 3-9 37 40 41 98; 7859) 120; 79709 10 24 71-74 77; 8074a 68; 8112i 1; 
8150i 1 126; 8165 10; 8249 b (8250) /8250e 161132 LYS tRNA: 8299k 126 

ATP 8: 8391e 31 120; 8515c27 ATP 6: 8592} 21 23; 8783g 6 7; 8852f (8854, 
8856) 109 110; 8994e60 61; 9009a 15 16; 9053f 52 91-93 95-102 112; 90701 5-7; 
9150j 7 103 104 CO lll: 9266e 28 49 67; 9294e 25; 9342e 27; 9380f 90; 9429 
56 57; 9553e 37 117; 9714e 98; 9746k 62; 97511 48; 9859g 86 GLY tRNA: 
10028a 126 ND3: 10066f 45; 100841 45; 10352a 110; 10364e 96; 10394c 1-29 
39-45 49 50 84-89 94-100 111 120-123 ARGtRNA: 10413a (10536) 16 46 56 
85 98 ND 4L: 10644k 19 20; 10689e 96; 10694a 104; 10725e 22 46 96 ND 4: 
10806g 1-7 50 83 124; 108931 34; 11146c 61 110; 171611 6 7; 11329e (11690) 2 
4, 17350a 15 37 110; 17806a 65; 11922] 7; 72026h 124 NDS5: 12345k 
(12350,12528) 88 127-134; 12406h 91-93; 12560a 19 52 93; 12795} 
(12798,12806,13374) 40 90 98 99; 12810k 5-7; 12925g 7; 12990a (13642) 27; 
13004) (13018,13182,13194) 20; 130319 1-6 9-16 18-23 26-35 37 39 40 42 43 
45-64 82 83 92 93; 13051e 23 38-41; 13068a 128; 13096k 26; 131001 107; 
13103g 7; 13208f 4 5 7 15 75 88; 132689 65; 13367b 40 88 90 99; 134041 8 9 
79; 136351 (13641) 39 115; 13702e 38-41; 14015a 14; 140501 (14366) 125 ND 6: 
14279e 1; 14279) 23; 14322a 92; 14385c 10; 14509a 68 123-126; 145671 68; 
14608c 5 CYT b: 14749e 87; 14869j 86 126; 149561 115; 15005g 56 57; 15172e 
10; 15195} (15221) 44; 15234g 113; 15238¢ 84 100; 15250c 14-18 20 21 24 25 
54 81 90 111; .75606a 65 127-134; 15723g 43; 15790j (16373) 82; 15883e 7 45 
58 59 112 THRtRNA: 15897k 97 129; 15907k 46; 159721 81; 15925) 15 81 90 
PRO tRNA: 15996c/16000g 48 ODLOOP: 16049k 15 22 46 74 81; 16065g 121 
122; 16089k 21; 16096k 119; 16125k 3 9 74 78 88-91 93 105 119 120 122; 
16778 26-29; 16208k 11 15 51 52 55; 162171 95; 162469 46; 16254a 127 132; 
16303k 36 58 59 82 89-93 111 112 131 132; 16310k 1-3 5 6 8-14 20 26-31 58 
73 78 82 94-97 99-102 120 125 126 133 134; 16389g9/16390b 11 35-41 82 89; 
16398e 5 11 34 37 66 118 119; 16490g 10; 16517e 2 11-13 15-24 26-29 51-55 
58 59 82 97 105-134. 


Fig. 3 a, Genealogical tree for 134 types of human mtDNA (133 restriction maps plus reference sequence); b, comprehensive list of polymorphic 
restriction sites used. The tree accounts for the site differences observed between restriction maps of these mtDNAs with 398 mutations. No 
other order of branching tested is more parsimonious than this one. This order of branching was obtained (using the computer program PAUP, 
designed by Dr David Swofford) by ignoring every site present in only one type of mtDNA or absent in only one type and confining attention 
to the remaining 93 polymorphic sites. The computer program produces an unrooted network, which we converted to a tree by placing the 
root (arrow) at the midpoint of the longest path connecting the two lineages (see ref. 58). The numbers refer to mtDNA types, no. 1 being 
from the aboriginal South African (!Kung) cell line (GM 3043), no. 45 being from the HeLa cell line and no. 110 being the published human 
sequence”’. Black bars, clusters of mtDNA types specific to a given geographic region; asterisks, mtDNA types found in more than 1 individual: 
type 134 was in six individuals, types 29, 65 and 80 each occurred thrice, and other types flagged with asterisks occurred twice. To place the 
nodes in the tree relative to the percent divergence scale, we took account of the differences observed at all 195 polymorphic sites. b, The 
numbering of sites is according to the published human sequence”*, with 12 restriction enzymes indicated by the following single letter code: 
a, Alul; b, Avall; c, Ddel; d, FnuDII; e, Haelll; f, Hhal; g, Hinfl; h, Hpal; i, Hpall; j, Mbol; k, Rsal; 1, Taql. Italicized sites are present 
in the indicated mtDNA types and nonitalicized sites are absent in the indicated types; parentheses refer to alternative placements of inferred 
sites; sites separated by a slash are polymorphic for two different restriction enzymes caused by a single inferred nucleotide substitution; 
letters in bold face refer to noncoding regions and genes for transfer RNA, ribosomal RNA and proteins (ND, NADH dehydrogenase; CO, 
cytochrome oxidase; ATP, adenosine triphosphatase; CYT b, cytochrome b). For example, 8j indicates a Mbol site beginning at nucleotide 
position 8 in the D loop that was not found in mtDNA type 1 but was present in type 2, etc. Note that since this site is not present in the 
reference sequence (type 110), the sequence beginning at position 8 is actually a semisite, differing from the MbolI recognition sequence at 
one position (see ref. 23 for a detailed description of the method of mapping such inferred sites). Not all sites were scored in all individuals: 
8j, 1484e and 7750c were not determined for types 1, 31, 59, 63, 68, and all of the New Guinea mtDNAs; mtDNAs 114 and 121 could not 
be typed with Rsal. The locations of some sites differ from those reported before”?*, as do the individuals in which some sites occur; these 
, revisions are based on re-examination of previously-studied mtDNAs. 


elsewhere*’. This rate is similar to previous estimates from 
animals as disparate as apes, monkeys, horses, rhinoceroses, 
mice, rats, birds and fishes'®. We therefore consider the above 
estimate of 2%-4% to be reasonable for humans, although 
additional comparative work is needed to obtain a more exact 
calibration. 

As Fig. 3 shows, the common ancestral mtDNA (type a) links 
mtDNA types that have diverged by an average of nearly 0.57%. 
Assuming a rate of 2%-4% per million years, this implies that 
the common ancestor of all surviving mtDNA types existed 
140,000-290,000 years ago. Similarly, ancestral types b-j may 
have existed 62,000-225,000 years ago (Table 3). 

When did the migrations from Africa take place? The oldest 
of the clusters of mtDNA types to contain no African members 
stems from ancestor c and included types 11-29 (Fig. 3). The 


apparent age of this cluster (calculated in Table 3) is 90,000- 
180,000 years. Its founders may have left Africa at about that 
time. However, it is equally possible that the exodus occurred 
as recently as 23-105 thousand years ago (Table 2). The mtDNA 
results cannot tell us exactly when these migrations took place. 


Other mtDNA studies 


Two previous studies of human mtDNA have included African 
individuals”’*; both support an African origin for the human 
mtDNA gene pool. Johnson et al”! surveyed ~40 restriction 
sites in each of 200 mtDNAs from Africa, Asia, Europe and the 
New World, and found 35 mtDNA types. This much smaller 
number of mtDNA types probably reflects the inability of their 
methods to distinguish between mtDNAs that differ by less than 
0.3% and may account for the greater clustering of mtDNA 
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% Sequence divergence 


Fig. 4 Genealogical tree relating the nucleotide sequences of D 
loops from seven human mtDNAs. This tree, which requires fewer 
mutations than any other branching order, was constructed by the 
PAUP computer program from the 900-bp sequences determined 
by Greenberg et al”® and was rooted by the midpoint method. 
Symbols: @, African origin (Black American); (, Caucasian. 


types by geographic origin that they observed. (By contrast, our 
methods distinguish between mtDNAs that differ by 0.03%.) 
Although Johnson et.al. favoured an Asian origin, they too 
found that Africans: possess. the greatest amount of mtDNA 
variability and that a midpoint rooting of their tree leads to an 
African origin. 

Greenberg et al2* sequenced the large noncoding region, 
which includes the displacement loop (D loop), from four 
Caucasians and three Black Americans. A parsimony tree for 
these seven D loop sequences, rooted by the midpoint method, 
appears in Fig. 4. This tree indicates (1) a high evolutionary rate 
for the D loop (at least five times faster than other other mtDNA 
regions), (2) a greater diversity among Black American D loop 
sequences, and (3) that the common ancestor was African. 


Nuclear DNA studies 


Estimates of genetic distance based on comparative studies of 
nuclear genes and their products differ in kind from mtDNA 
estimates. The latter are based on the actual number of muta- 
tional differences between mtDNA genomes, while the former 
rely on differences in the frequencies of molecular variants 
measured. between and within populations. Gene frequencies 
can be influenced by recombination, genetic drift, selection, and 
migration, so the direct relationship found between time and 
mutational distance for mtDNA would not be expected for 
genetic distances based on nuclear DNA. But. studies based on 
polymorphic blood groups, red cell enzymes, and serum proteins 
show that (1) differences between racial groups are-smaller than 
those within such groups and (2) the largest gene frequency 
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differences are between Africans and other populations, suggest- 
ing an African origin for the human nuclear gene pool'®'*>. 
More recent studies of restriction site polymorphisms in nuclear 
DNA"**¢-” support these conclusions. 


Relation to fossil record 


Our tentative interpretation of the tree (Fig. 3) and the associated 
time scale (Table 3)-fits with one view of the fossil record: that 
the transformation of archaic to anatomically modern forms of 
Homo sapiens’ occurred first in Africa*’*, about 100,000- 
140,000 years ago, and that all present-day humans are descen- 
dants of that African population. Archaeologists have observed 
that blades were in common use in Africa 80-90 thousand years 
ago, long before they replaced flake tools in Asia or Europe***’. 
But the agreement between our molecular view and the evidence 
from palacoanthropology and archaeology should be treated 
cautiously for two reasons. First, there is much uncertainty about 
the ages of these remains. Second, our placement of the common 
ancestor of all human mtDNA diversity in Africa 140,000- 
280,000: years ago need not imply that the transformation to 
anatomically modern Homo sapiens occurred in Africa at this 
time. The mtDNA data tell us nothing of the contributions to 
this transformation by the genetic and cultural traits of males 
and females whose mtDNA became extinct. 

An alternative view of human evolution rests on evidence that 
Homo has been present in Asia as well as in Africa for at least 
one million years*® and holds that the transformation of archaic 
to anatomically modern humans occurred in parallel in different 
parts of the Old World® 49 “This hypothesis leads us to expect 
genetic differences of great antiquity within widely separated 
parts of the modern pool of mtDNAs. It is hard to reconcile 
the mtDNA results with this hypothesis. The greatest divergences 
within clusters specific to non-African parts of the World corre- 
spond to times of only 90,000-180,000 years. This might imply 
that the early Asian Homo (such as Java man and Peking man) 
contributed no surviving mtDNA lineages to the gene pool of 
our species. Consistent with this implication are features, found 
recently in the skeletons of the ancient Asian forms, that make 
it unlikely that Asian erectus was ancestral to Homo sapiens °°, 
Perhaps the non-African erectus population was replaced by 
‘sapiens migrants from Africa; incomplete fossils indicating the 
possible presence of early modern humans in western Asia at 
Zuttiyeh (75,000-150,000 years ago) and Qafzeh (50, 000-70, 000 
years.ago) might reflect these first migrations** 

If there was hybridization between the resident archaic forms 
in Asia and anatomically modern forms emerging from Africa, 
we should expect to find extremely divergent types of mtDNA 
in present-day Asians, more divergent than any mtDNA found . 
in Africa. There is no evidence for these types of mt DNA among 
the. Asians studied?'**-*°, Although such archaic types of 
mtDNA could have been lost from the hybridizing population, 
the probability of mtDNA lineages becoming ‘extinct in an 
expanding population is low’’. Thus we propose that Homo 





Table 3 Ancestors, lineages and extents of divergence in the geneaological tree for 134 types of human mtDNA 





No. of descendant lineages or clusters specific to a region 


Ancestor Total Africa Asia Australia Europe’ N. Guinea % divergence Age* 
a 7 1 0 0 0. 0 0.57 143-285 
b 2 0 1 0 0 0 0.45 112-225 
c 20 0 7 3 1 3 0.43 108-215 
d 2 0 0 1 1 0 0.39 98-195 
e 14 - 2 2 4 2 0 0.34 85-170 
f 19 1 7 4 4 1 0.30 75-150 
g 10 2 3 2 2 1 0.28 70-140 
h 30 2 4. 0 15 1 0.27 68-135 
i 8 1 0 0 6 0 0.26 65-130 
j 22 1 3 1 5 1 0.25 62-125 

All 134 10 27 15 36 7 — _ 





* Assuming that the mtDNA divergence rate is 2-4% per million years’ 
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erectus in Asia was replaced without much mixing with the 
invading Homo sapiens from Africa. 


Conclusions and prospects 


Studies of mtDNA suggest a view of how, where and when 
modern humans arose that fits with one interpretation of 
evidence from ancient-human bones and tools. More extensive 
molecular comparisons are needed to improve our rooting of 
the mtDNA tree and the calibration of the rate of mtDNA 
divergence within the human species. This may provide a more 
reliable time scale for the spread of human populations and 
better estmates of the number of maternal lineages involved in 
founding the non-African populations. 

It is also important to’ obtain more quantitative estimates of 
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Apparent faster-than-light motion has been detected with very long 
baseline interferometry (VLBI) in the central regions of ~12 
extragalactic radio sources”, These are typically strong, core- 
dominated objects and, so far the only weak-cored classical double 
radio,source to show the effect was 3C179 (ref. 3). Unified source 
models*, which invoke relativistic motion ata small angle to the 
observer’s line of sight, predict that large superluminal: motion 
should be rarely seen in a randomly oriented sample of such 
objects. These models have motivated statistical studies of quasar 
samples’°, designed. to detect superluminal objects and, sub- 
sequently, ‘to compare the distribution of velocities with model 
predictions. Here we report the detection of superluminal motion 
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the overall extent of nuclear DNA diversity in both human and 
African ape populations. By comparing the nuclear and 
mitochondrial DNA diversities, it may be possible to find out 
whether a transient or prolonged bottleneck in population size 
accompanied the origin of our species!*. ‘Then a fuller inter- 
action between palaeoanthropology, archaeology and molecu- 
lar biology will allow a deeper analysis of how our species 
arose. 
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in the quasar 3C263. The milliarc second structure of this source 
consists of two nearly unresolved components, whose separation 
is increasing at a rate of 0.06+ 0.02 mare s yr ". This corresponds 
to an expansion speed of 2.70.9 c (Hy =S55kms~ * Mpe™" 1%= 
0.05). This quasar is the weakest superluminal source found so 
far, and there are indications. that superluminal: motion occurs 
frequently in this class of object. - 

The quasar 3C263 (=1;137+660; m,=16, variable’, z= 
0.652'') was earlier mapped with arc second resolution at 
5 GHz", at 2.7 and 8.1 GHz'”, and at 1.42 GHz’. It has a steep 
overall spectrum™ and consists of a weak core (~130 mJy at 
5 GHz) with inverted spectrum and two dominant steep-spec- 
trum lobes with separation of ~45” in position angle (p.a.) 110°. 


’ There is no evidence for radio variability of the total source flux 


density at 5GHz'> or at 10.7 GHz"®, nor is variability of the 
central component detected at 15 GHz (J. F. C. Wardle, personal 
communication). Earlier VLBI observations of the .coré at 
5 GHz" and at 10.7 GHz’ showed an extended core structure 
aligned with the outer double lobes. 

We observed the core of 3C263 at thrée epochs with a trans- 
atlantic VLB array operating at a frequency of 10.7 GHz (left- 
hand circular polarization). These observations were made. at 
intervals of ~1 yr, each time with up to five antennas {Table 1). 
The Mark III recording system'® was used in mode A or B (56 
or 28 MHz bandwidth). Data were recorded for ~12 hin scans 


